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<~i . Abstract 

In this article, I present an elementary introduction to the theory 
of gravitational waves. This article is meant for students who have 
had an exposure to general relativity, but, results from general rela- 
^ I tivity used in the main discussion has been derived and discussed in 

the appendices. The weak gravitational field approximation is first 
^.^1 considered and the linearized Einstein's equations are obtained. We 

• discuss the plane wave solutions to these equations and consider the 

Q-i| transverse-traceless (TT) gauge. We then discuss the motion of test 

particles in the presence of a gravitational wave and their polarization. 
The method of Green's functions is applied to obtain the solutions to 
the linearized field equations in presence of a nonrclativistic, isolated 
""V . source. 
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1 Introduction 



In the past few years, research on the detection of gravitational waves has 
assumed new directions. Efforts are now underway to detect gravitational 
radiation from astrophysical sources, thereby enabling researchers to pos- 
sess an additional tool to study the universe (See Q for a recent review). 
According to Newton's theory of gravitation, the binary period of two point 
masses (e.g., two stars) moving in a bound orbit is strictly a constant quan- 
tity. However, Einstein's general theory of relativity predicts that two stars 
revolving around each other in a bound orbit suffer accelerations, and, as a 
result, gravitational radiation is generated. Gravitational waves carry away 
energy and momentum at the expense of the orbital decay of two stars, 
thereby causing the stars to gradually spiral towards each other and giving 
rise to shorter and shorter periods. This anticipated decrease of the orbital 
period of a binary pulsar was first observed in PSR 1913+16 by Taylor and 
Weisberg {^). The observation supported the idea of gravitational radia- 
tion first propounded in 1916 by Einstein in the Proceedings of the Royal 
Prussian Academy of Knowledge. Einstein showed that the first order con- 
tributon to the gravitational radiation must be quadrupolar in a particular 
coordinate system. Two years later, he extended his theory to all coordinate 
systems. 

The weak nature of gravitational radiation makes it very difficult to 
design a sensitive detector. Filtering out the noisy background to isolate 
the useful signal requires great technical expertise, itself a field of research. 
Various gravitational wave detectors are fully /partially operational and we 
expect a certain result to appear from the observations in the near future. 

This article gives an elementary introduction to the theory of gravita- 
tional waves. Important topics in general relativity including a brief intro- 
duction to tensors and a derivation of Einstein's field equations are discussed 
in the appendices. We first discuss the weak gravitational field approxima- 
tion and obtain the linearized Einstein's field equations. We then discuss 
the plane wave solutions to these equations in vacuum and the restriction on 
them due to the transverse-traceless (TT) gauge. The motion of particles in 
the presence of gravitational waves and their polarization is then discussed 
in brief. We conclude by applying the method of Green's functions to show 
that gravitational radiation from matter at large distances is predominantly 
quadrupolar in nature. 
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2 The weak gravitational field approximation 



Einstein's theory of general relativity leads to Newtonian gravity in the limit 
when the gravitational field is weak & static and the particles in the gravita- 
tional field move slowly. We now consider a less restrictive situation where 
the gravitational field is weak but not static, and there are no restrictions 
on the motion of particles in the gravitational field. In the absence of grav- 
ity, space-time is flat and is characterised by the Minkowski metric, r]^^. A 
weak gravitational field can be considered as a small 'perturbation' on the 
flat Minkowski metric[^, 

gfiu = rjp.u + h^^, < 1 (1) 

Such coordinate systems are often called Lorentz coordinate systems. Indices 
of any tensor can be raised or lowered using r]^'^ or r]^i, respectively as 
the corrections would be of higher order in the perturbation, h^y. We can 
therefore write, 

gi"^ = r/^^ - h^"' (2) 

Under a background Lorentz transformation (Qj), the perturbation trans- 
forms second-rank tensor: 

/i„;3 = A„''A/3" V (3) 

The equations obeyed by the perturbation, /i^j^, are obtained by writing the 
Einstein's equations to first order. To the first order, the affine connection 
(See App. A) is, 

= lv^n^^.K, + d.h^p - V] + (^{h^) (4) 
Therefore, the Riemann curvature tensor reduces to 

Rfiupa = Vfj.xdpT — ri^xdfj^ (5) 
The Ricci tensor is obtained to the first order as 



jl ^(1) = 1 r 



(6) 



where, □ = r]^^d\dp is the D'Alembertian in flat space-time. Contracting 
again with r/^'^, the Ricci scalar is obtained as 

R = dxd^h^'' - ah (7) 
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The Einstein tensor, G^u, in the hmit of weak gravitational field is 

(8) 

The linearised Einstein field equations are then 

G^, = STTGrf"" (9) 

We can't expect the field equations (9) to have unique solutions as any solu- 
tion to these equations will not remain invariant under a 'gauge' transforma- 
tion. As a result, equations (9) will have infinitely many solutions. In other 
words, the decomposition (1) of Qf^i, in the weak gravitational field approx- 
imation does not completely specify the coordinate system in space-time. 
When we have a system that is invariant under a gauge transformation, we 
fix the gauge and work in a selected coordinate system. One such coordinate 
system is the harmonic coordinate system (|5|). The gauge condition is 

5^"r\, = o (10) 

In the weak field limit, this condition reduces to 

dxh\ = \d,h (11) 

This condition is often called the Lorentz gauge. In this selected gauge, the 
linearized Einstein equations simplify to, 

□V - Ivi^u^h = -WnGT"" (12) 

The 'trace-reversed' perturbation, /i^^,, is defined as ([^), 

hf,u = hf,u - ^ili^,yh (13) 
The harmonic gauge condition further reduces to 

d^h\ = (14) 

The Einstein equations are then 

□ /i^^ = -IGTrGT'^'' (15) 
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3 Plane- wave solutions and the transverse-traceless 
(TT) gauge 

From the field equations in the weak-field limit, eqns.(15), we obtain the 
linearised field equations in vacuum, 

□V = (16) 

The vacuum equations for hf^i^ are similar to the wave equations in electro- 
magnetism. These equations admit the plane-wave solutions, 

h-nu = Ai^i,exp{LkaX°') (17) 

where, Aj^^, is a constant, symmetric, rank-2 tensor and ka is a constant 
four-vector known as the wave vector. Plugging in the solution (17) into the 
equation (16), we obtain the condition 

kak"" = (18) 

This implies that equation (17) gives a solution to the wave equation (16) 
if ka is null; that is, tangent to the world line of a photon. This shows 
that gravitational waves propagate at the speed of light. The time-like 
component of the wave vector is often referred to as the frequency of the 
wave. The four- vector, k^ is usually written as k^ = (a;, k). Since ka is null, 
it means that, 

= \kf (19) 

This is often referred to as the dispersion relation for the gravitational wave. 
We can specify the plane wave with a number of independent parameters; 
10 from the coefficients, A^i, and three from the null vector, A;^. Using the 
harmonic gauge condition (14), we obtain, 

kaA""^ = (20) 

This imposes a restriction on : it is orthogonal (transverse) to ka- The 
number of independent components of is thus reduced to six. We have 
to impose a gauge condition too as any coordinate transformation of the 
form 

x"' = a;" + ^°(x^) (21) 
will leave the harmonic coordinate condition 

□x'^ = (22) 
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satisfied as long as 



Let us choose a solution 



(23) 



(24) 



to the wave equation (23) for Cq, are constant coefficients. We claim 
that this remaining freedom allows us to convert from the old constants, 



.(old) 



(new) 



ifiu , to a set of new constants, A^u 

-- 



A (new) 
^,1 



such that 
(traceless) 



(25) 



and 



(26) 

where, is some fixed four-velocity, that is, any constant time-like unit 
vector we wish to choose. The equations (20), (25) and (26) together are 
called the transverse traceless (TT) gauge conditions ([^). Thus, we have 
used up all the available gauge freedom and the remaining components of 
A^u must be physically important. The trace condition (25) implies that 



h 



TT 

a/3 



h 



TT 



(27) 



Let us now consider a background Lorentz transformation in which the 
vector C/" is the time basis vector = Then eqn.(26) implies that 
A^o = for all /i. Let us orient the coordinate axes so that the wave 
is travelling along the z-direction, — > (a;,0, 0,u;). Then with eqn.(26), 
eqn.(20) implies that A^z = for all a. Thus, A^ in matrix form is 



A 



TT 
a/3 



( ° 








M 





A 

■^xx 


A 

^xy 








A 

■'^xy 


_A 
■'^xx 





I 








0/ 



(28) 



4 Polarization of gravitational waves 

In this section, we consider the effect of gravitational waves on free particles. 
Consider some particles described by a single velocity field, U" and a sepa- 
ration vector, (^°. Then, the separation vector obeys the geodesic equation 
(See App. A) 



(29) 
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where, = dx^ /dr is the four- velocity of the two particles. We consider the 
lowest-order (flat-space) components of only since any corrections to 
that depend on h^i, will give rise to terms second order in the perturbation in 
the above equation. Therefore, = (1, 0, 0, 0) and initially C,^ = (0, e, 0, 0). 
Then to first order in h^y, eqn. (29) reduces to 

Using the definition of the Riemann tensor, we can show that in the TT 

gauge, 



nx 

^ 0x0 — 


RxOxO = 


_^uTT 
2 '''xx,00 


0x0 — 


RyOxO = 


1 TT 


OyO — 


RyOyO = 


2"'TO,00 



-R\xO (31) 

All other independent components vanish. This means that two particles 
initially separated in the x-direction have a separation vector which obeys 
the equation 

Similarly, two particles initially separated by e in the y-direction obey the 
equations 

g2^y _ 1 52 _ 1 52 

52c- _1 d\TT 

5*2 -2^dt^ ^y ^ ' 

We can now use these equations to describe the polarization of a gravita- 
tional wave. Let us consider a ring of particles initially at rest as in Fig. 
1(a). Suppose a wave with /i^J / Oj^xJ — ^ ^^^^ them. The particles 
respond to the wave as shown in Fig. 1(b). First the particles along the 
x-dircction come towards each other and then move away from each other 
as reverses sign. This is often called -|- polarization. If the wave had 
^xy 7^ ^1 but, /i^J = hyy = 0, then the particles respond as shown in Fig. 
1(c). This is known as x polarization. Since h'^y and /i^J are independent, 
the figures 1(b) and 1(c) demonstrate the existence of two different states 
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of polarisation. The two states of polarisation are oriented at an angle of 
45" to each other unlike in electromagnetic waves were the two states of 
polarization. 



5 Generation of gravitational waves 

In section III, we obtained the plane wave solutions to the linearized Ein- 
stein's equations in vacuum, eqns.(16). To obtain the solution of the lin- 
earised equations (15), we will use the Green's function method. The Green's 
function, G{x^ — yf^), of the D'Alembertian operator □, is the solution of 
the wave equation in the presence of a delta function source: 

□ Gixf" - yf") = (x^ - y") (34) 

where 5^^^ is the four-dimensional Dirac delta function. The general solution 
to the linearized Einstein's equations (15) can be written using the Green's 
function as 

V(a;") = -167rG j d'y G(x" - y")r^.(y") (35) 

The solutions to the eqn.(34) are called advanced or retarded according as 

they represent waves travelling backward or forward in time, respectively. 
We are interested in the retarded Green's function as it represents the net 
effect of signals from the past of the point under consideration. It is given 

by 

G{x^^ - y^) = -^-^—^6 [|x - y| - (a;" - y")] x ^(xO - y^) (36) 

where, x = {x^,x^,x'^) and y = {y^,y'^,y^) and |x — y| = [dij{x'^ — y^){x^ — 
y')]^/^. 9{x^ — y^y is the Heaviside unit step function, it equals 1 when 
x^ > y^, and equals otherwise. Using the relation (36) in (35), we can 
perform the integral over y^ with the help of the delta function, 



V(i,x) =4G / d^y ^^r^,(t- |x-y|,y) (37) 
where t = The quantity 

tR = t - |x - y| (38) 
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is called the retarded time. From the expression (37) for hn^,, we observe that 
the disturbance in the gravitational field at (t, x) is a sum of the influences 
from the energy and momentum sources at the point (tR,y) on the past 
light cone. 

Using the expression (37), we now consider the gravitational radiation 
emitted by an isolated far away source consisting of very slowly moving 
particles (the spatial extent of the source is negligible compared to the dis- 
tance between the source and the observer). The Fourier transform of the 
perturbation h^i, is 

^^^(a;,x) = — ^ f dt exp(ia;t) h^„{t,y:) (39) 
y/27r J 

Using the expression (37) for hni, {t, x) , we get 

V = d^y exp(ta;|x - y|) (40) 



Under the assumption that the spatial extent of the source is negligible com- 
pared to the distance between the source and the observer, we can replace 
the term exp(ta;|x — y|)/|x — y| in (40) by exp(ia;R)/R. Therefore, 

^^.(c.,x)=4G^^P|^^ jd^yf,,{u,y) (41) 
The harmonic gauge condition (14) in Fourier space is 

- f ^ ap 

d^h '"'{t,x) = df, dco h exp(-ia;i) = (42) 
Separating out the space and time components, 

do dio h {LO,x)exp{—iLot) + di dio h (w, x)exp(— twt) = (43) 

Or, 

LLoh = dih (44) 

Thus, in eqn.(41), we need to consider the space-like components of hni,{u, y). 
Consider, 

j d^y dk {y%j) = J d^y (dky') f^^ + j d^y y' {d^f^^) 
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On using Gauss' theorem, we obtain, 

I d^y r^icj, y) = - / d^y y' (^^r^O ^^^^ 

Consider the Fourier space version of the conservation equation for T^^'^, viz., 
d^T^^^{t,yi) = 0. Separating the time and space components of the Fourier 
transform equation, we have, 

diT" = LuT^"" (46) 

Therefore, 

j d?y r^iiv, y)=iivj cfy y^ f^^ = j d'y (y^ f^^ + ^ f^') (47) 
Consider 



d'y di (y, y,f^') = j d% [(diy^) y^ + (dw^) f] f^'+ J d'y y^ y^ (s^f o') 

(48) 

Simphfying the equation above, we obtain for the left hand side 

I d^y (y* ro^- +y^fo*) + j d'y y' y^ {^dif^') 

Since the term on the left hand side of eqn.(47) is a surface term, it vanishes 
and we obtain 

J d^y (y'f^^ + y^rO») =- j d^y y' y^ (dif^') (49) 
Using the equations (46) and (48), we can write, 

j d^y r^(a;,y) = '-^jd'y di (yVf^') (50) 
Using the eqn(45), we can write, 

I d^y r\u,y) = - y/ d^'y yW (51) 

We define the quadrupole moment tensor of the energy density of the source 
as 

QijiLo) = ^jd'y yY f^'^u;, y) (52) 
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In terms of the quadrupole moment tensor, we have 



/ 



d^yf\uj,y) = -'^qij{uj) (53) 



Therefore, the solution (41) becomes 



i / \ exp(6a;R) / ^'^ ~ , A ,rA\ 

hij (u;, x) = 4G I - —Qij (u) 1 (54) 

Simphfying further, 

= 2 Gui'^ 

hij{u,x) = - - — — exp(6a;R) qij{uj) (55) 
o K 

Taking the Fourier transform of eqn.(54), and simphfying, we finaUy obtain 
for the perturbation 

2G d 

where, iR, = t — |x — y| is the retarded time. In the expression (54), we see 
that the gravitational wave produced by an isolated, monochromatic and 
non-relativistic source is therefore proportional to the second derivative of 
the quadrupole moment of the energy density at the point where the past 
light cone of the observer intersects the cone. The quadrupolar nature of 
the wave shows itself by the production of shear in the particle distribution, 
and there is zero average translational motion. The leading contribution to 
electromagnetic radiation comes from the changing dipole moment of the 
charge density. This remarkable difference in the nature of gravitational 
and electromagnetic radiation arises from the fact that the centre of mass 
of an isolated system can't oscillate freely but the centre of charge of a 
charge distribution can. The quadrupole momnet of a system is generally 
smaller than the dipole moment and hence gravitational waves are weaker 
than electromagnetic waves. 



6 Epilogue 

This lecture note on gravitational waves leaves several topics untouched. 
There are a number of good references on gravitation where the inquisitive 
reader can find more about gravitational waves and their detection. I have 
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freely drawn from various sources and I don't claim any originality in this 
work. I hope I have been able to incite some interest in the reader about a 
topic on which there is a dearth of literature. 
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Appendix A: Some topics in general theory of relativity 



An event in relativity is characterised by a set of coordinates {t, x, y, z) 
in a definite coordinate system. Transformations between the coordinates 
of an event observed in two different reference frames are called Lorentz 
transformations. These transformations mix up space and time and hence 
the coordinates arc redefined so that all of them have dimensions of length. 
We write x^ = ct, = x, 
vector {x^ , x^ , x'^ , x^ 



y,x^ = z and a general component of a four 
as x^. A Lorentz transformation is then written as 



x"^ 



(57) 



where, 



A = 



7 
-Pi 





-/37 

7 








1/ 



(58) 



At this point, it is useful to note the Einstein summation convention: when- 
ever an index appears as a subscript and as a superscript in an expression, 
we sum over all values taken by the index. Under a Lorentz transformation, 
the spacetime interval — (ct)^ + + + remains invariant. The length 
of a four-vector is given by 



2 + (^2)2 ^ (^3)2 



(59) 



We never extract a square root of the expression (59) since |x| can be nega- 
tive. Four-vectors that have negative length are called time-like, while those 
with positive lengths are called space-like. Four-vectors with zero length are 
called null. The notion of "norm" of a four-vector is introduced with the 
help of the Minkowski metric: 



/ -1 \ 

10 

10 

\ 1 / 



(60) 



Then, we have. 



X^tJui^X 



(61) 
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There are two kinds of vectors that are classified in the way they transform 
under a Lorentz transformation: 



Contravariant = ^x^ 

Covariant = A^^Xi, (62) 

Vectors are tensors of rank one. r]'^'^ (rjni,) is called the metric tensor; it 
is a tensor of rank two. There are other higher rank tensors which we 
will encounter later. If two coordinate systems are linked by a Lorentz 
transformation as: 

x"' = A'' ^xi" (63) 

then, multiplying both sides of the equation above by A^^ and differentiat- 
ing, we get, 

A. ^ (64) 



dx' ^ 

Therefore, we see that 



' (65) 



dx' ^' ^ dx 
Thus, 

^l d d d d \ 

transforms as a covariant vector. The differential operates on tensors to 
yield higher-rank tensors. A scalar s can be constructed using the Minkowski 
metric and two four-vectors u^^ and v"^ as: 

s = Vi^^u^v" (67) 

A scalar is an invariant quantity under Lorentz transformations. Using the 
chain rule, 

dx' " = — — dx"" (68) 



we have. 



If we define 



dx'' 



^^'"dx"' dx'^ ) 



u"'v'^ (69) 



then. 



dx^ dx 

5«A = r/,.^^ (70) 
s = g^x u'^v'^ (71) 
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g^x is called the metric tensor; it is a symmetric, second-rank tensor. 

To follow the motion of a freely falling particle, an inertial coordinate 
system is sufficient. In an inertial frame, a particle at rest will remain 
so if no forces act on it. There is a frame where particles move with a 
uniform velocity. This is the frame which falls freely in a gravitational field. 
Since this frame accelerates at the same rate as the free particles do, it 
follows that all such particles will maintain a uniform velocity relative to this 
frame. Uniform gravitational fields are equivalent to frames that accelerate 
uniformly relative to inertial frames. This is the Principle of Equivalence 
between gravity and acceleration. The principle just stated is known as the 
Weak Equivalence Principle because it only refers to gravity. 

In treating the motion of a particle in the presence of gravity, we define 
the Christoffel symbol or the affine connection as 

1 — y ^ « -I- „ „ yi-i) 



r plays the same role for the gravitational field as the field strength tensor 
does for the electromagnetic field. Using the definition of affine connection, 
we can obtain the expression for the covariant derivative of a tensor: 

^«^'^ = ^ + r\.^" (73) 

It is straightforward to conclude that the covariant derivative of the metric 
tensor vanishes. The concept of "parallel transport" of a vector has im- 
portant implications. We can't define globally parallel vector fields. We 
can define local parallelism. In the Euclidean space, a straight line is the 
only curve that parallel transports its own tangent vector. In curved space, 
we can draw "nearly" straight lines by demanding parallel transport of the 
tangent vector. These "lines" are called geodesies. A geodesic is a curve of 
extremal length between any two points. The equation of a geodesic is 

(Px'^ ^„ dx^ dxf^ 

The parameter A is called an affine parameter. A curve having the same path 
as a geodesic but parametrised by a non-affinc parameter is not a geodesic 
curve. The Riemannian curvature tensor is defined as 
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In a "flat" space, 

= (76) 

Geodesies in a flat space maintain their separation; those in curved spaces 
don't. The equation obeyed by the separation vector in a vector field V 
is 

vv vv c = R^^a.y^ y" (77) 

If wc differentiate the Ricmannian curvature tensor and permute the indices, 
we obtain the Bianchi identity: 

dxRajBiiv + 9t,Rai3\ij. + dfj^Ra/Bi^x = (78) 

Since in an inertial coordinate system the affine connection vanishes, the 
equation above is equivalent to one with partials replaced by covariant 
derivatives. The Ricci tensor is defined as 

Rap = = (79) 

It is a symmetric second rank tensor. The Ricci scalar (also known as scalar 
curvature) is obtained by a further contraction, 

R = Rf^i^ (80) 

The stress- energy tensor (also called energy-momentum tensor) is defined as 
the flux of the a-momentum across a surface of constant x^. In component 

form, we have: 

1. r'^^ = Energy density = p 

2. = Energy flux (Energy may be transmitted by heat cinduction) 

3. T'^ = Momentum density (Even if the particles don't carry momen- 
tum, if heat is being conducted, then the energy will carry momentum) 

4. T'-' = Momentum flux (Also called stress) 
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Appendix B: The Einstein field equation 



The curvature of space-time is necessary and sufficient to describe grav- 
ity. The latter can be shown by considering the Newtonian hmit of the 
geodesic equation. We require that 

• the particles are moving slowly with respect to the speed of light, 

• the gravitational field is weak so that it may be considered as a per- 
turbation of flat space, and, 

• the gravitational field is static. 

In this limit, the geodesic equation changes to, 

^^ + rS.(|)-o ,81) 

The affine connection also simplifies to 

r^o = -^ff^'5Affoo (82) 

In the weak gravitational field limit, we can lower or raise the indices of a 
tensor using the Minkowskian flat metric, e.g., 

v'^'Kp = h% (83) 
Then, the affine connection is written as 

C = -\n^^dxhoo (84) 
The geodesic equation then reduces to 

The space components of the above equation are 

(fx' 1 , dt 



Or, 



7;i-j-rdihoo (86) 



J2 i 1 
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The concept of an inertial mass arises from Newton's second law: 



f = m^a (88) 

According to the the law of gravitation, the gravitational force exerted on an 
abject is proportional to the gradient of a scalar field ^, known as the scalar 
gravitational potential. The constant of proportionality is the gravitational 
mass, nig: 

fg = -rUgV^ (89) 

According to the Weak Equivalence Principle, the inertial and gravitational 
masses are the same, 

mi = mg (90) 

And, hence, 

a = -V$ (91) 

Comparing equations (86) and (91), we find that they are the same if we 
identify, 

/loo = -2$ (92) 

Thus, 

500 = -(1 + 2$) (93) 

The curvature of space-time is sufficient to describe gravity in the Newto- 
nian limit as along as the metric takes the form (93). All the basic laws of 
Physics, beyond those governing freely-falling particles adapt to the curva- 
ture of space-time (that is, to the presence of gravity) when we are work- 
ing in Riemannian normal coordinates. The tensorial form of any law is 
coordinate-independent and hence, translating a law into the language of 
tensors (that is, to replace the partial derivatives by the covariant deriva- 
tives), we will have an universal law which holds in all coordinate systems. 
This procedure is sometimes called the Principle of Equivalence. For exam- 
ple, the conservation equation for the energy-momentum tensor T>^^ in fiat 
space-time, viz., 

is immediately adapted to the curved space-time as, 

D^T^''' = (95) 

This equation expresses the conservation of energy in the presence of a 
gravitational field. We can now introduce Einstein's field equations which 
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governs how the metric responds to energy and momentum. We would hke 
to derive an equation which will supercede the Poisson equation for the 
Newtonian potential: 

= -AttGp (96) 

where, = S^^didj is the Laplacian in space and p is the mass density. 
A relativistic generalisation of this equation must have a tensorial form so 
that the law is valid in all coordinate systems. The tensorial counterpart of 
the mass density is the energy-momentum tensor, T^*^. The gravitational 
potential should get replaced by the metric. Thus, we guess that our new 
equation will have T^^'^ set proportional to some tensor which is second-order 
in the derivatives of the metric, 

T'^" = (97) 

where, Af^i^ is the tensor to be found. The requirements on the equation 
above are:- 

• By definition, the R.H.S must be a second-rank tensor. 

• It must contain terms linear in the second derivatives or quadratic in 
the first derivative of the metric. 

• The R.H.S must be symmetric in /i and v as T^" is symmetric. 

• Since T'*'^ is conserved, the R.H.S must also be conserved. 
The first two conditions require the right hand side to be of the form 

aR^^ + (3Rg^^ = T^^ (98) 

where R^^ is the Ricci tensor, R is the scalar curvature and a &: /3 are 
constants. This choice is symmetric in jjL and v and hence satisfies the third 
condition. Prom the last condition, we obtain 

g'"'D„{aR^, + l3Rg^,) = Q (99) 

This equation can't be satisfied for arbitrary values of a and 13. This equa- 
tion holds only if a/ (5 is fixed. As a consequence of the Bianchi identity, 
viz., 

Di'R^,^ = ^D^R (100) 
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we choose, 

/3 = -ia (101) 
With this choice, the equation (42) becomes 

aiRf,, - ^Rgf,,) = T^"" (102) 

In the weak field hmit, 

ffoo ~ -2$ (103) 
the 00-component of the equation(42), viz., 

- aV^c/oo = Too ^ 2aV2$ = p (104) 
Compare this result with Newtons equation (40), we obtain. 

Thus, we obtain the Einstein field equations in their final form as 

Ri^u - IrQi^u = SnGT'^'' (106) 
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Figures 




Fig 1(a) 



Fi gure 1 ; The initial configuration of test particles on a circle before a gravitational wave hits 
them. 




t 

Figure 2: Displacement of test particles caused by the passage of a gravitational wave with the 
+ polcirization. The two states are separated by a phase difference of n. 
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Fig. 1(c) 



Figure 3: Displacement of test particles caused by the passage of a gravitational wave with the 
X polarization. 
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